arXiv:1508.04894vl [astro-ph.SR] 20 Aug 2015 


The cluster ages experiment (CASE). VII. Analysis of two 
eclipsing binaries in the globular cluster NGC 6362 

J. Kaluzny^t, I. B. Thompson^, A. Dotter^, M. Rozyczka^, A. Schwarzenberg-Czerny^ 

^Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland; 

mnr@camk. edu.pl, alex@camk. edu.pl 

“^Observatories of the Carnegie Institution of Washington, 813 Santa Barbara Street, Pasadena, 

CA 91101-1292, USA; ian@obs.carnegiescience.edu 

^Research School of Astronomy and Astrophysics, Australian National University Canberra, 

Australia; aaron.dotter@gmail.com 


August 21, 2015 


ABSTRACT 

We use photometric and spectroscopic observations of the detached eclips¬ 
ing binaries V40 and V41 in the globular cluster NGC 6362 to derive masses, 
radii, and luminosities of the component stars. The orbital periods of these 
systems are 5.30 and 17.89 d, respectively. The measured masses of the pri¬ 
mary and secondary components (Mp, M^) are (0.8337±0.0063, 0.7947T0.0048) 
M0 for V40 and (0.8215T0.0058, 0.7280T0.0047) M© for V41. The measured 
radii (Rp, R,) are (1.3253±0.0075, 0.997±0.013) Rq for V40 and (1.0739±0.0048, 
0.7307±0.0046) Rq for V41. Based on the derived luminosities, we hnd that the 
distance modulus of the cluster is 14.74±0.04 mag - in good agreement with 
14.72 mag obtained from CMD htting. We compare the absolute parameters of 
component stars with theoretical isochrones in mass-radius and mass-luminosity 
diagrams. For assumed abundances [Fe/H] = -1.07, [a/Fe] = 0.4, and Y = 0.25 
we hnd the most probable age of V40 to be 11.7±0.2 Gyr, compatible with the 
age of the cluster derived from CMD htting (12.5±0.5 Gyr). V41 seems to be 
markedly younger than V40. If independently conhrmed, this result will suggest 
that V41 belongs to the younger of the two stellar populations recently discovered 
in NGC 6362. The orbits of both systems are eccentric. Given the orbital pe¬ 
riod and age of V40, its orbit should have been tidally circularized some Gyr 
ago. The observed eccentricity is most likely the result of a relatively recent close 
stellar encounter. 
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1. Introduction 


The Gluster AgeS Experiment (GASE) is a long term project conducted mainly at 
Las Gampanas Observatory with a principal aim of the detection and detailed study of 
detached eclipsing binaries (DEBs) in nearby globular clusters (GGs). The main goal is 
to measure masses, luminosities and radii of DEB components with a precision of better 
than on e per cent in order to determine GG ages and distances, and to test stellar evolution 
m odels flKahiznv et al.N2nn5h. The methods and assumptions we employ follow the ideas 
of iPaczvhskil fjl997() and [Thompson et al.l (1200 ll) . Thus far, we have presented results for 
seven DEBs from four clusters: u Gen, 47 Tuc, M4 and M55. These were the hrst and, 
to the best of our knowledge, the only direct measurements of the fundamental parameters 
of main-sequence and subgian t stars in GGs. More d etails on the project and references to 
earlier papers can be found in iKaluznv et al.l (j2014al) . 


The present paper is devoted to an analysis of two DEBs, V40 and V41, belonging 
to the globular cluster NGG 6362. Radial velocity and photometric observations are used 
to derive masses, luminosities, and radii of the component stars. In Section [2] we describe 
the photometric observations and determine orbital ephemerides. Section E] presents spec¬ 
troscopic observations and radial-velocity measurements. The combined photometric and 
spectroscopic solutions for the orbital elements and the resulting component parameters are 
obtained in Section |H In Section |5l we compare the derived parameters to Dartmouth stel¬ 
lar evolution models, with an emphasis on estimating the age of the binaries. Finally, in 
Section E] we summarize our Endings. 


fDeceased 

*This paper includes data gathered with the 6.5-in Magellan Baade and Clay Telescopes, and the 2.5-ni 
du Pont Telescope located at Las Campanas Observatory, Chile. 

























3 


2. Photometric observations 

Photometric observations analyzed in this paper were obtained at Las Campanas Ob¬ 
servatory between 1995 and 2011 using the 2.5-m du Pont telescope equipped with the TEK5 
and SITE2K CCD cameras, and the 1.0-m Swope telescope equipped with the SITES CCD 
camera. The scale was 0.26 arcsec/pixel and 0.44 arcsec/pixel for the du Pont and Swope 
telescopes, respectively. For each telescope the same set oi BV Liters was used for the entire 
duration of the project. The first likely eclipse events of V40 and V41 were observed in 
June 1991. For the next eight years observations were conducted in a survey mode with the 
aim of detecting more eclipses and to establish epher nerides for bo t h vari ables. The periodic 
variability of V40 and V41 was hrmly established bv iMazur et al.l (1l999l) . Starting with the 
2001 season, we attempted to cover predicted eclipses of V40 or V41. Several eclipses were 
observed with both telescopes. The light curve analysis in Section 0] is based exclusively 
on the du Pont data. The less accurate Swope photometry was used solely to improve the 
ephernerides (the difference in accuracy between the two telescopes resulted from the differ¬ 
ence in the light-gathering power and the poorer sampling of the point spread function of 
the Swope telescope). 

The core and half-lig ht radii of NGC 6362 are = 70" and = 123", respectively 
flHarrislll996. 2010 editionl) . V40 is located in the core region at a distance of d = 59" from 
the cluster center. V41 resides at the half-light radius with a distance of d = 121". The 
du Pont images show no evidence of unresolved visual companions to either binary. In the 
case of V40 this is conhrmed by the examination of HST/ACS images obtained within the 


6362 ( 

Zloczewski et al. 

2012 : 

Kaluznv et al. 

201 4bh 

and V41 can be found in 

Mazur et al. 

(1999 

)■ 


Coordinates and Ending charts for V40 


The light curves presented h ere were measured with profile fitting photometry extracted 
with the Daophot/Allstar codes flStetsonlll98h. Il990l) . and transformed to the standard John¬ 
son system based on observations of Landolt f 19921) standards. A more detailed description 


of the data and the reduction procedures can be found in a paper pr esenting the photometry 
of 69 variables from the field of NGC 6362 f Kaluznv et al. 2014b). We note that, in gen¬ 


eral, more accurate photometry of crowded fields can be obtained with an image subtraction 
technique. However, profile fitting photometry is less vulnerable to systematic errors, and 
as such it is a better option whenever the star is free from blending. 
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2.1. Ephemerides 


Several eclipse events of V40 were covered densely enongh to determine individnal times 
of minima. They were fonnd using t he KWALEX code wh i ch is based on an improved 
version of the algorithm presented by Kwee fc van Woerden f IQsd) (see Appendix). The 
linear ephemerides: 


^ 2453915.50966(25) + 5.2961725(14) x E 
= 2453891.82853(14) + 5.29617486(99) x E (1) 


proved to be an adequate £t to the moments of minima listed in Table [U The orbit of V40 
is eccentric. Since the secondary (shallower) minimum was observed more times than the 
primary one, we decided to choose it as a base of the phase count. Thus, we set = 0.5, 

which necessarily results in ^ 0. 


Altogether, 13 eclipses of V41 were observed between August 2007 and September 2011. 
However, ingress and egress were only covered for a few of them. This prevented us from a 
classical period study based on the determination of individual moments of minima and the 
0-C technique. To find the ephemeris, a combined V light curve including data from all sea- 
sons and both tele scopes was used. This was fitted with the help of the JKTEBOP v34 cod^ 
(1Southworthll2013l and references therein), allowing only for variations of P and moment of 
the primary eclipse Tq. Mode 8 of the code was used, which enabled a robust determination 
of the errors of the two parameters with the help of Monte Carlo simulations. The following 
linear ephemeris 


^ 2454263.65242(25) + 17.8888441(40) x E (2) 


was obtained. The light curves phased with ephemerides given by equations ([T]) and ([2]), and 
used for the subsequent analysis, are shown in Fig. [1] 


3. Spectroscopic observations 


The spectra were taken with the MIKE echelle spectrograph flBernstein et al.ll2003l) on 
the Magellan Baade and Clay 6.5-m telescopes, using a 0.7" slit which provided a resolution 
R ~40,000. A typical observation consisted of two 1800-second exposures of the target, 
flanking an exposure of a thorium-argon hollow-cathode lamp. The raw spectra were re¬ 
duced with the pipeline software written by Dan Kelson, following the approach outlined 


^described in detail at and available from littp://www.astro.keele.ac.uk/jkt/codes/jktebop.html 
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in 


KelsonI (120031) . For the further processing, the IRAF package ECHELL E was used. The 


veloc ities were measured using software based on the TODCOR algorithm fjZucker &: Mazeh 


19941) . kindly made availabl e by G uillermo Torres. Synthetic echelle-resolution spectra from 


the library of ICoelho et al.l (120061) were used as velocity templates. These templates were 
interpolated to the values of log^f and Te// derived from the photometric solution (see Sec¬ 
tion Hj) and Gaussian-smoothed to match the resolution of the observed spectra. For the 
interpolation in chemical composition we assumed [Fe/H] = —l-O? and fa/Fe] = 0.4, i.e. 
values adopted for the cluster as a whole (ICaretta et al.ll2009l: IPotter et al.N2010[) . These 
values are not well determined (see the discussion in Sect. II, however the results of the 
velocity measurements are insensitive to minor changes in [Fe/H] or [a/Fej. 


Radial velocities were measured on the wavelength intervals 4120-4320 A and 4350-4600 
A, covering the region of the the best signal-to-noise ratio while avoiding the Hy line. The 
hnal velocities, obtained by taking the mean of these two measurements, are given in Tables 
[2] a nd [3l which list the heliocentric Julian dates at mid-exposure, velocities of the primary 
and secondary components, and the orbital phases of the observations calculated according 
to the ephemerides given by Equations (1) and (2). An estimate of the error in the velocity 
measurements is the rms of half of the differences between velocities obtained from the two 
wavelength intervals. The corresponding values are 0.63, 0.59, 0.23 and 0.52 km/s for the 
primary and secondary of V40 and the primary and secondary of V41, respectively. They 
are consistent with the rms values of the residuals from orbital hts given in Table lH 


4 . Analysis of light and velocity curves 


We assume that reddening is uniform for NGC 6362. iHarrisI (11996. 2010 edition! ) lists 
E{B — V) = 0.09 as measured by West and Zinn (1984) based on the Q39 index. The map 
of total galactic extinction by Schlegel et al. (1998) predicts E{B — V) = 0.075 mag for the 
position of NGG 6362. This value agrees with an est imate resulting from isochrone htting 
applied to HST/AGS photometry of the cluster bv iDotter et ahl (120101) . They obtained 
E{Q — 8) = 0.070, equivalent to E(B — V) = 0.071. This value is also com patible with 


values obtained by iKovacs fc Walkerl (I2OOII) {E{B — V)= 0.073 ± 0.007) and lOlech et al. 


(I2OOII) {E{B — V) = 0.08 ± 0.01) from an analysis of RR Lyr stars in NGG 6 362. We adopt 
E(B — V) = 0.075 and neglect the diff erential reddening, w hich according to (IBonatto et al. 


20131 1 is 0.025 mag on the average. (IBonatto et al.l 120131) comment that this differential 
reddening may result from uncorrected zero-point variations in the original photometry. 


We employed a quadratic formula for limb darkening, and calculated coefficients it¬ 
eratively together with effective temperatures of the components of both systems. The 
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temperatures were deriv ed from dereddened B — V colors using the empirical calibra tion of 


Casaerande et al. 

(2010 

). The coefficients were then interpolated from the tables of 

Claret 


total V and B magnitudes at maximum light and light ratios derived from the light curve 
solution. We note that a slightly more recent calibration by Sousa et al. (2011) results 
in temperatures higher by 8 K on average, however the agreement between calibrations is 
excellent over the relevant range of colors. 


We checked that using a square root or logarithmic formula for the limb darkening 
instead of a quadratic formula had a negligible impact on the photometric solution. In 
particular, relative radii derived from the three two-parameter formulas differed by less than 
~ 0 . 2 cr), whereas using a linear limb darkening formula resulted in a slightly poorer £t, and 
caused the relative radii to differ by ~lcr from those obtained with the quadratic formula. 
We conclude that, for the stars analyzed here, two-parameter formulate provide a better 
description of the limb darkening than the linear law. 


Light and radial velocity curves of V40 and V41 were analyzed with the JKTEBOP v34 
code, which allows for a simultaneous solution of one light curve and two radial velocity 
curves (one for the primary and one for the secondary component). The Id-curve was solved 
for the following parameters: inclination i, sum of fractional radii Vp and r^, ratio of the radii 
Ts/rp] surface brightness ratio S (secondary to primary), luminosity ratio Lg/Lp, ecos(a;) and 
esin(a;) where e is the orbital eccentricity and oj the longitude of the periastron, amplitudes 
of radial velocity curves Kp and Kg, and systemic velocity 7 . Errors of stellar parameters 
were derived using a Monte Carlo approach implemented in mode 8 of the JKTLD code. We 
performed 10000 Monte Carlo simulations for each of V40 and V41. 


Since the 5-data were of poorer quality than those collected in the Id-band, we decided 
to use them solely to a measurement of the luminosity ratio {Lg/Lp)B- This was accomplished 
by fixing all parameters but the surface brightness ratio Sb- However, the formal error of 
Sb calculated by JKTEBOP would have been seriously underestimated if uncertainties of 
other parameters had not been taken into acco unt. To avoid this problem we used the 
PHOEBE package vO.Sla fjPrsa &: Zwitteiil2005l) . which allows for a simultaneous solution 
of several light curves but can be unstable when attempting to solve light and velocity 
curves simultaneously. We initialized the package with the photometric solution found by 
JKTEBOP, loaded 5 Id curves, and, allowing all geometrical parameters to vary, found the 
standard error of using another Monte Carlo procedure written in PHOEBE-scripter. 
The procedure replaces the observed light curves Bq and Vq with the htted values Bf and 


^Written by John Southworth and available at www.astro.keele.ac.uk/jkt/codes/jktld.html 
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Vf , generates Gaussian perturbations 6Bf and 6Vf such that the standard deviation of 
the perturbation is equal to the standard deviation of the corresponding residuals shown in 
Figure [2], and performs PHOEBE iterations on Bf + 6Bf and Vf + 6Vf. As before, 10000 
Monte Carlo simulations were performed for each binary. In our opinion such a hybrid 
approach based on both the JKTEBOP and PHOEBE codes is optimal when analyzing 
data sets similar to those presented in this paper. 


The quality of the final fits is illustrated in Figures [2] and [21 Orbital and photometric 
parameters derived for both systems are presented in Table 01 and absolute parameters 
of the components are listed in Table [5l The errors of the effective temperatures include 
uncertainties of observed colors, reddening, and photometric zero points (0.015 mag for 
V and B — V). The bolometric luminosities are calculated from radii and temperatures 
adopting Lq* = 3.827 x 10^^ ergsj^ To derive absolute visual mag n itude s, we used bolometric 
corrections based on atmosphere models of ICastelli fc Kuruczl (120041) which ranged from 
—0.07 mag for the primary of V40 to —0.10 mag for the secondary of V41. The absolute 
visual magnitudes along with the H-magnitudes derived from light curve solutions allowed 
for a determination of the apparent distance moduli listed in Table O The location of the 
components of V40 and V41 on the CMD of NGC 6362 is shown in Figure 01 All four stars 
are placed on the main sequence: three of them at the turnoff region, and the fourth one 
(the secondary of V41) well below the turnoff. 


5. Age and distance analysis 


The fundamental parameters obtained from the analysis of V40 and V41 allow us to 
derive the ages of the components of these systems, using theoretical isochrones. However, 
before a model comparison is made it is necessary to review the information available on the 
chemical composition of NGG 6362 since the model-based ages are sensitive to the adopted 
values of helium abundance, [Fe/H] and [a/Fe]. 


Unfortunately, the metallicity of NGG 6362 is not well known. A literature search re- 
ve als that determinatio n s of Fe/H] from high resolution spectra are limited to the study 
by Garretta fc Gratton f l997 h who found [Fe/H] = —0.96 and [Fe/H] = —0.97 for two 


red giants. iGaretta et ah! (j2009h used this determination along with some older estimates 
based on the integrated photometry of NGG 6362 to derive [Fe/H] = —1.07 ± 0 .05 on 
their newly propose d metallicity scale. For the sake of completeness, we note that [Harris 
f[l996. 2010 editioru ) listed [Fe/H] = —0.99 - a value obtained by combining that o nGaretta et ah 


^ https://sites.google.com/site/mamajeksstarnotes/basic-astronomical-data-for -the-sun 
























020091) with [Fe/H] = —0.74 ± 0.05 derived by Ge isler et al. (1997) based on Washington 


photometry of 10 cluster giants. lOlech et al.l (120011) ob t ained a helium content Y=0.292 from 


the analysis of RRc stars. Finally, iDalessandro et al.l (120141) presented convincing eyidence 
of multiple populations in NGC 6362. We adopt a hducial composition of [Fe/H]=—1.07, 
[a/Fe]=+0.4 and Y=0.25, but will investigate the dependence on these parameters when 
deriving the ages of the component stars. 


In a study of three binary s ystems in M4 flKaluzny et al.l 1201311 we used two sets of the- 
oretical isochrones: Dart mouth flDotter et al.l 120081 henceforth DSED) and Victoria-Regina 
(IVandenBerg et al.ll2012[) . Since the differences in the derived ages turned out to be smaller 
than the uncertainties imposed by the observational errors, the present analysis is based 
on DSED isochrones only. Figure 0 shows DSED isochrones for ages ranging from 9-Gyr 
to 13-Gyr plotted on mass-radius and mass-luminosity diagrams, together with values for 
V40 and V41. Both diagrams indicate that V41 may be younger than V40, and the mass- 
radius diagram suggests that the secondary of V41 is younger than the primary (or, if both 
components are of the same age, the secondary is too small compared to the primary). 

A detailed comparison of the ages of the individual components, which requires a quanti¬ 
tative age est imate from the data a nd the isochrones, was undertaken following the approach 
developed bv iKaluznv et al.l (120131) to study the binary systems in M4. Briefly, the (M — R) 
or (M — L) plane is densely sampled in both dimensions for each component separately, and 
the age of an isochrone matching the properties of the sampled point is found. The age 
ensembles generated in this way are visualized in the form of mass-age diagrams (MADs). 
The MAD of a given component is composed of two ellipses derived from (M — R) and 
(M — L) relations. Each ellipse dehnes a range of ages achieved by stars whose parameters 
differ from the corresponding entries (M, R) or (M, L) in Table l5]by at most Er or El, where 
— (t^M/MqY -1- a{R/RQY for the (M — R) plane and likewise for the (M — L) plane; 
a^M/Mo), a^R/Ro) and a^L/Lo) being the standard errors of M, R and L from Tabled 
The MADs of the components of V40 and V41 are shown in Figure [61 

The age of each component was calculated as the weighted mean of the ages derived 
from the two corresponding MADs. For V40 we obtained 11.7±0.3 Gyr for the primary 
and 11.5±0.4 Gyr for the secondary. The difference is 0.2T0.5 Gyr. For V41 the ages 
are 10.7±0.3 Gyr and 9.2T0.6 Gyr, respectively, for the primary and the secondary. The 
difference is 1.5±0.7 Gyr. Possible origins of this two-sigma effect are discussed in Section |6l 
We estimated the age of each system from the weighted mean of the ages derived from each 
pair of correspo nding MADs. For V 40 we obtained 11.7±0.2 Gyr, in agreement with 12.5±0.5 
Gyr derived by Potter et al. ( 2010 ) from GMD htting. For V41 we obtained 10.4T0.3 Gyr, 
which suggests that this system is signihcantly younger than V40, with a formal difference 
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of 1.3±0.4 Gyr. 


To independently estimate the age of the cluster as a whole, we htted the CMD of 
NGC 6362 with DSED isochrones (see Figure [7]). Since uncertainties of stellar models grow 
larger with evolutionary time, when fitting we assigned the largest weight to main sequence 
and turnoff stars. This results in an age estimate of 12.5±0.75 Gyr - again in agreement 
with the result of iDotter et al.l (120101) . which is based on an independent set of photomet¬ 
ric data. The GMD fit y ielded a distance modulus of 14.72 mag - larger than 14. 55 mag 
derived bv iDotter et al.l (120101) . but in good agreement with 14.68 mag listed by iHarris 
( 1996. 2010 edition ). The reddening derived from the GMD-fit {E{B — V) = 0.061 mag) is 
slightly lower than the value of 0.075 mag adopted in Sect. HI but decreasing [Fe/H] from 
the fiducial -1.07 to -1.15 brings it back to the adopted value (we checked that such a small 
change in metallicity does not influence the parameters listed in Tables H] and E] in any sig¬ 
nificant way). The distance moduli of all four stars forming V40 and V41 agree with each 
other (see Table [5]), and the mean weighted modulus of 14.74±0.04 mag is consistent with 
that obtained from GMD fitting. 

To sum up, the age derived from the MADs of V40 components is compatible with that 
derived from GMD htting, however the age of the primary of V41 is 1.0±0.4 Gyr younger than 
the age of V40, with the secondary of V41 being another 1.5±0.7 Gyr younger. The problem 
looks much milder in the (M — L) diagram, in which all the four components are marginally 
coeval within l-cx limits. We stress, however, that the [M — L) diagram and associated MADs 
must be given a smaller weight beca use the luminositie s in Table E] are sensitive to adopted 
chemical composition (via Teff, see IDotter et al.l 120101) and possible differential reddening, 
whereas the radii are entirely reddening and composition-independent. 


6. Discussion and summary 


We derived absolute parameters of the components of V40 and V41, two detached 
eclipsing binaries located near the main sequence turnoff of NGG 6362. The accuracy of 
mass and radii measurements is generally better than 1 per cent, with the only exception 
being the radius of the secondary in V40 which is determined with an accuracy of 1.3 per 
cent. The luminosities of the components w ere found using effective temperatures estimated 
from {B — V) - calibrations compiled bv iGasagrande et al.l (120101) . While their errors are 
rather large (between 7 and 10 per cent), all the four distance moduli have nearly the same 
value, and the formally calculated mean modulus of 14.74T0.04 mag agrees very well with 
14.72 mag derived from the comparison of DSED isochrones with the GMD of the cluster. 
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Based on the parameters of V40 and DSED isochrones for Y=0.25, [Fe/H] =—1.07 and 
[a/Fe] =+0.4, we derived an age of 11.7+0.2 Gyr. The good agreement of this resnlt with the 
age of NGC 6362 derived from isochrone comparison with the GMD of the clus ter is enconr- 


aging. This is particularly the case given the age of the cluster as obtained by iDotter et ah 


fl2010l ) based on an independent photometry (12.5+0.5 Gyr). However, one should keep in 
mind that the uncer tainty of 0.2 Gyr is a purely statistical one, and, as pointed out by 
Kaluznv et al.l (1201311 . it should be increased by ~0.85 Gyr to account for the sensitivity of 
stellar age to uncertainties in the chemical composition. 

The situation with V41 is less clear. Isochrone hts to (M — R) and (M — L) diagrams 
of this system suggest that it is signihcantly (by 1.3 + 0.4 Gyr) younger than V40. However, 
given the recent discovery of two stellar populations in NGG 6362, such a possibility does not 
seem to be entirely imp lausible. Upon examining GMDs obtained from HST photometry, 
Dalessandro et al.l (1201411 found that the red giant branch of the cluster splits into two sub¬ 
branches separated by about 0.1 mag in {U — V), while such effect is not observed in (V — I). 
Their Figure 1 explains the split by a difference in light-element content. While they use 
isochrones of the same age, which is reasonable given the relatively weak dependence of RGB 
location on time, they clearly identify the Na-poor populations with the hrst generation of 
stars, and the Na-rich population with the second one ’’formed during the hrst few ~100 
Myr of the cluster life from an intra-cluster medium polluted by the hrst generation”. We 
observed similar ehects (i.e. a large split in [U — V) and a much smaller one in (U —/)) using 
same-age DSED isochrones dihering by ~0.1 dex in [Fe/H], or ~0.3 dex in [a/Fe], or ~0.1 in 
Y. The RGB also splits when the populations have the same composition but diher in age by 
4-5 Gyr, however an age diherence this large causes the subgiant branch to split even more 
widely th an the RGB, while t he op posite is observed in NGG 6362. Thus, we corroborate the 
results of lDalessandro et al.l (120141 1 and their scenario involving two star formation episodes, 
with the second-generation stars slightly enriched in metals (and possibly in helium). 

The discrepancy of the ages of V41 components (1.5 + 0.7 Gyr) may be regarded as 
statistically insignihcant. However we decided to study the possible origin of the difference by 
varying chemical composition parameters Y, [Fe/H] and [a/Fe] (only one of them was varied 
at a time), and looking for isochrones joining locations of the primary and the secondary 
in the (M — R) and (M — L) diagrams. We found acceptable (albeit only marginally) hts 
with Y = 0.26 at an age of 9.5 Gyr; [a/Fe] = 0.2 at an age of 9.7 Gyr, and [Fe/H] = -0.90 
at an age of 11.3 Gyr. This suggests that V41 may have formed from material enriched in 
products of stellar nucleosynthesis, and thus it may indeed be younger than V40. We note 
parenthetically that the results of our htting experiments are in line with the Endings of 
Kaluznv et al.l (120131) . who, while discussing (M — R) diagrams for binaries in M4, wrote 
“increasing Y decreases the age derived from the (M — R) plane while increasing [Fe/H] or 
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[a/Fe] increases the age.” 

V40 is also somewhat peculiar because of the eccentricity of its orbit. At an age of ~12 
Gyr and a period of 5.3 d the orbit shou ld long ago have been fully circularized by tidal 
friction flMazehl 120081: iMathieu et ahl 120041) . There is no evidence for the presence of a third 
body in this system, which means that it mus t have undergone a close encounter during the 
last few Gyr. According to iDalessandro et al.l (120141) NGG 6362 may have had a complicated 
dynamical history, losing perhaps as much as 80% of its original mass. Since close encounters 
account for a few per cent of the total mass loss (M. Giersz, private communication), the 
very existence of the nonzero eccentricity of V40 lends some support to this conjecture. 

While V40 and V41 seem to belong to distinct populations, we feel it is too early to 
claim a direct detection of non-coevality among stars in globular clusters. It is clear, however, 
that DEBs offer the potential to test theories describing the origin of multiple populations in 
GGs, and to verify estimates of age differences between populations based on GMD-fitting. 
For the latter, both the theoretical pred ictions and fitting pr ocedures yield values ranging 
from ~ 10^ yr to several hundred Myr (INardiello et ahl 120151 and references therein), and 
any independent constraints are highly desirable. The most obvious test for the reality of the 
age difference between our two systems would be to d etermine the ch emical composition of 
their components using disentangling software (see e.g. lHadravall2009l) . The existing spectra 
do not have adequate quality to measure abundances, but, given the brightnesses of the 
components and the orbital periods, it should be possible to improve the S/N ratio. 

There is a l so roo m for improving the accuracy of parameter measurements. As noted by 
Kaluznv et al.l (120131) . the errors in the masses of the components originate almost entirely 
from the orbital solution, which may only be improved by taking additional spectra (prefer¬ 
ably with the same instrument). Admittedly, this would require a substantial observational 
effort, as doubling of the present se t of radial velocity data would lead to an improvement 
of only 33% in the mass estimates (iThompson et al.l 120101) . The errors in the radii of the 
components are dominated in turn by the photometric solution, whose quality depends on 
the accuracy of the photometry. We estimate that an accuracy of 0.002 - 0.003 mag in V 
would reduce the errors of the radii by 50%. This goal could be easily achieved on a 6-8 m 
class telescope equipped with a camera capable of good PSF sampling. Finally, IR photom¬ 
etry should be obtained to reduce errors in luminosity with the help of existing accurate 
calibrations linking {V — K) color to R-band surface brightness. 


This series of papers is dedicated to the memory of Bohdan Paczyhski. JK and MR were 
supported by the grant DEG-2012/05/B/ST9/03931 from the Polish National Science Gen- 
ter. AD received support from the Australian Research Gouncil under grant FL110100012. 
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Appendix: Kwalex: A method for timing eclipses 


Kwee fc van WoerdenI fll956l) proposed an analyis of the timing of eclipses by fitting 
them with symmetric cnrves, so that any asymmetries wonld be reflected in observed-minns- 
calcnlated {O — C) deviations. Their classical method does have its drawbacks. As far as 
the algorithm is concerned: (i) the interpolation of a light curve into a somewhat arbitrary 
evenly spaced time mesh yields a slight dependence of a hnal estimate of the central time 
on its initial approximation and (ii) the resulting estimate of the eclipse centre depends on 
the initial selection of data points and on the initial guess of the eclipse centre. There are 
also statistical drawbacks: (hi) interpolation effectively assigns uneven weights to points and 
(iv) the effective model, a polyline, suffers from excess number of parameters describing all 
nodes. Such a waste of degrees of freedom is known in statistics to decrease accuracy. 

Let us adopt the zero point of the tirn e scal e as the current estimate of the eclipse 
centre Ti"'\ Following Kwee fc van Woerden ( 1956 ) we £t both the input light curve and its 
reflexion with respect to Tc'^\ The novelty here is that we employ the orthogonal polynomial 
approximation of the combined light curve (e.g. iRalstonI 1196511 .rather than on a poly-line 
interpolation. In doing so we optimize the sensitivity of our method by keeping the number 
of model paramete rs at a statist ically justihable minimum (Occam razor). By virtue of the 
Fisher lemma (e.g. lBrandtlll9751) . parameters of the orthogonal model are uncorrelated with 
noise hence they yield the minimum variance model estimate. 


Next, again following iKwee fc van WoerdenI fjl956l) . we freeze the model curve and shift 
by ±5 and calculate of O — O' for the time shifts —5, 0, -l-(5. These 3 values of are 
htted with a parabola 

y^(t) = at^ + bt + c (3) 


and the location of its minimum = —b/2a may be adopted as the next approximation 

of However, to assure convergence, we trim time steps dTc = —Tc"'^ and care¬ 

fully select 6 to avoid going too far away from the solution, yet possibly ensuring bracketing 
it within ±5. Such orthogonal polynomials retain only even powers and are optimal both in 
terms of degrees of freedom and numerical procedure relying on orthogonal projection. By 
the incorporation of weight coefficients in the scalar product our algorithm is able to account 
for uneven errors of data. 


Complications, both conceptual and algorithmic, arise due to our ultimate aim of mini¬ 
mizing the dependence of the hnal result on somewhat accidental selection of the initial data 
points. Discrete selection/modification of observations may result in numerical instability 
and/or non-converging limit cycles. To minimize such danger we resort to smooth down¬ 
weighting of points near edge of eclipse by a symmetric Fermi-Dirac distribution function 
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(e.g. 


Kitte]|ll996[) : 


fit) = 


e s +1 


(4) 


The final weight of data measured at time ti with an error Uj is Wi = f{ti)/a‘f. As our time 
zero-point moved during iterations towards the final solution, so our weights vary between 
iteration s. We roughly followed princ iples of the minimization algorithm by simulated an¬ 
nealing fjOtten fc van GinnekenI Il989l ). At the start of iterations we adopt a rather crude 
model curve for m = 2 and smooth weighting 6 = A/10. Note, that in statistical physics 6 
would be proportional to the temperature. Only after successful progress of our algorithm do 
we tighten our data sample by decreasing 6 (i.e. ’lower temperature’) and refine our model 
by increasing m. Also, we trimmed the time step dTc < A/10. 


After each iteration we calculate two Tc error estimates: S and a. The error S cor¬ 
responds to such a model shift S = T, that grows over its minimum value Xmin by its 
variance \/M, where d = n — m/2 — 1 is the number of degrees of freedom. The error a is 
defined similarly, for growth by 1. For this purpose has to be renormalized by a factor 
C such, that E{Cx^} = d. For an approximate normalization we substitute the minimum 
value of the parabola in Eq. ([3]) for its expectation, so that C = (4ac — b‘^)/Aad. S and a 
are calculated by solving of Eq. (I3l). In this way we obtain: = \/^Xmin/(^ similarly 

for cr^: 




^/2 


4ac — 6^ 
4a^\/d 


4ac — 6^ 


4a‘^d 


(5) 

( 6 ) 


To force a positive result we apply Eq. ([3]) to log rather then to y^. In statistical terms, S 
corresponds to the maximum reasonable uncertainty of Tc, were other parameters strongly 
correlated with it and a is its likely error, assuming weak correlation. 


During iterations we keep 6 as large as 5 = 2S. Thus 6 decreases for an improved model 
fit at a moderate rate. The polynomial order is raised in steps of 2, up to its maximum 
allowed value, only after decrease of corrections of Tc to dTc = < O.OIS. We 

terminate iterations when m reaches its selected maximum value and time steps become 
small dTc < eA, where e is our tolerance parameter. Our result is the final eclipse centre 
Tc and its likely error a. Not e, that for correlated observations errors need to be increased 
fjSchwarzenberg-Czernyl Il99ll ) . 
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Table 1. 


Times of minima for V40 


E 

HJD-2450000 

a 

0-C 

-480.5 

1349.66400 

0.00129 

0.00060 

-274.5 

2440.67681 

0.00083 

-0.00019 

0.5 

3891.82809 

0.00027 

0.00045 

3.5 

3907.71733 

0.00024 

-0.00027 

16.5 

3976.56802 

0.00139 

-0.00069 

77.5 

4299.63424 

0.00026 

-0.00024 

138.5 

4622.70052 

0.00025 

0.00015 

276.5 

5353.57276 

0.00031 

0.00003 

0 . 

3915.50987 

0.00028 

-0.00022 

1 . 

3920.80387 

0.00096 

0.00196 

139.0 

4651.67750 

0.00035 

0.00014 

335.0 

5689.72753 

0.00039 

-0.00007 
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Table 2. Velocity observations of V40 


HJD-2450000 

Vp [km s 

Vs [km s 

phase 

2870.5433 

52.73 

-78.06 

0.665 

3183.6775 

57.18 

-85.40 

0.790 

3201.5720 

-74.19 

52.16 

0.169 

3206.5845 

-65.72 

41.55 

0.115 

3210.5848 

35.22 

-64.64 

0.871 

3517.7456 

36.05 

-63.77 

0.867 

3521.7538 

40.02 

-66.70 

0.624 

3580.5437 

60.61 

-87.52 

0.725 

3585.6468 

56.59 

-84.07 

0.688 

3586.6453 

33.28 

-59.01 

0.877 

3815.8955 

-74.88 

52.41 

0.163 

3816.8897 

-62.77 

41.06 

0.350 

3877.8517 

38.86 

-66.45 

0.861 

3890.7002 

-74.76 

55.72 

0.287 

3935.6890 

57.56 

-86.79 

0.782 

3937.6503 

-72.90 

50.94 

0.152 

3938.5635 

-68.71 

47.47 

0.324 

3991.5267 

-68.84 

47.33 

0.325 

4258.7016 

59.64 

-88.13 

0.771 

4258.7465 

57.69 

-86.47 

0.780 

4258.7904 

58.16 

-85.71 

0.788 

4329.5454 

-72.97 

49.67 

0.148 

4648.6241 

-49.91 

27.55 

0.395 

4656.6819 

15.95 

-41.39 

0.916 
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Table 3. Velocity observations of V41 


HJD-2450000 

Vp [km s 

Vs [km s 

phase 

2871.5571 

-32.73 

8.93 

0.181 

2872.5519 

-38.27 

15.99 

0.236 

3179.7609 

-47.66 

27.81 

0.410 

3206.5336 

20.39 

-52.18 

0.906 

3517.7940 

-43.64 

23.06 

0.306 

3520.7452 

-46.33 

26.13 

0.471 

3521.7074 

-40.89 

20.15 

0.525 

3580.6885 

42.56 

-73.79 

0.822 

3581.6731 

30.56 

-59.28 

0.877 

3582.5946 

16.30 

-44.91 

0.928 

3891.6858 

-34.69 

14.52 

0.207 

3937.6050 

44.58 

-76.40 

0.774 

3938.6590 

40.59 

-73.04 

0.833 

3990.5131 

35.50 

-66.03 

0.731 

4329.6373 

15.02 

-42.95 

0.689 

4648.7184 

-41.08 

19.93 

0.525 

4649.6742 

-30.79 

8.04 

0.579 

4649.7192 

-30.27 

7.02 

0.581 

4652.7191 

40.41 

-73.29 

0.749 

4655.7155 

19.30 

-48.17 

0.917 

4967.7970 

-47.13 

25.65 

0.362 

5012.6510 

30.80 

-62.41 

0.870 

5038.5921 

-44.61 

23.99 

0.320 

6842.6401 

-30.67 

7.75 

0.167 

6844.6397 

-41.48 

20.30 

0.279 

6845.6556 

-46.30 

24.06 

0.336 
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Table 4. Orbital and photometric parameters®" of the systems V40 and V41 


Parameter V40 V41 


7 (km s“^) 
Kp (km s“^) 
Kg (km s“^) 
ap (km s“^) 
as (km s“^) 

^(Ro) 

e 

uj (deg) 
i (deg) 


Sb 

i^Lp/Ls')y 
{^Lp/L s') B 

o-rms(R) (mmag) 
o-rms(R) (mmag) 
Vp (mag)'’ 

Vs (mag)*® 

{B - V)p (mag)’ 
{B - V)s (mag)’ 


-12.34(11) 

70.18(17) 

73.63(25) 

0.72 

1.07 

15.037(32) 
0.05054(73) 
206.6(17) 
88.22(7) 
0.08817(45) 
0.06627(87) 
0.9815(35) 
0.9686(18) 
1.803(54) 
1.8272(57 ) 
11 
13 

18.698(13)(20) 

19.338(22)(27) 

0.542(18)(23) 

0.556(31)(34) 


-12.55(8) 

46.59(14) 

52.57(17) 

0.53 

0.65 

33.294(71) 

0.3125(13) 

321.23(33) 

89.547(20) 

0.03225(13) 

0.02195(12) 

0.7333(63) 

0.669(11) 

2.964(27) 

3.264(44) 

12 

16 

19.089(7)(17) 

20.274(10)(18) 

0.550(10)(18) 

0.650(16)(22) 


^Numbers in parentheses are the errors of the last 
signihcant digit (s). 

’’For V and {B — V) both the internal error (from 
the photometric solution and profile photometry) and 
the total error is given, the latter including 0.015 mag 
uncertainty of the zero point of the magnitude scale. 
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Absolute parameters^ of the components of V40 and V41 


Parameter 

V40 

V41 

Mp (Mo) 

0.8337(63) 

0.8215(58) 

M, (Mo) 

0.7947(48) 

0.7280(47) 

Rp (R-©) 

1.3253(77) 

1.0739(48) 

Rs (R©) 

0.997(13) 

0.7307(46) 

Tr(K) 

6156(123) 

6124(121) 

r. (K) 

6100(158) 

5747(122) 

ir (Le) 

2.27(18) 

1.46(11) 

(Ls) 

1.24(13) 

0.524(45) 

Mvp (mag) 

3.92(8) 

4.40(8) 

Mvs (mag) 

4.59(11) 

5.54(9) 

(m - M)vp (mag) 

14.78(8) 

14.69(8) 

(m - M)vs (mag) 

14.75(9) 

14.73(8) 

Agcp (Gyr) 

11.7(3) 

10.7(3) 

Ages (Gyr) 

11.5(4) 

9.2(6) 

Agcb (Gyr) 

11.7(2) 

10.4(3) 


^Numbers in parentheses are the errors of 
the last signihcant digit (s). The age of the 
binary Agcb is a weighted mean of component 
ages. 
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19.4 

0.48 0.5 0.52 

Phase 

Fig. 1.— Light curves adopted for the analysis: V40 (left column) and V41 (right column). 
In each column the complete curve is shown in the top panel, followed by zoomed primary 
minimum and secondary minimum. The 5-curves are shifted by 0.3 mag for clarity. 
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0.48 0.5 0.52 0.65 0.66 

Phase Phase 


Fig. 2.— Residuals from the fits to the light curves of V40 (left) and V41 (right). The 
panels are arranged in the same way as the two lower rows of Fig. [H In each panel, the 
lower sequence represents the R-residuals, and the upper one the R-residual offset by 0.1 
mag for clarity. 
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Fig. 3.— Observed velocities of the components of V40 (left) and V41 (right) shown together 
with the hnal hts. Filled and open symbols denote primary and secondary components, 
respectively (note that for V40 phase 0.0 corresponds to the secondary minimum). Errors of 
individual points are smaller than the symbols themselves. Lower panels display residuals 
from the hts. 
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Fig. 4.— The color-magnitude diagram of NGC 6362 with l ocations of V41, y42 ari d their 
components. Data for the background stars are taken from IZloczewski et ahl (120121 ) (only 
proper motion members of the cluster are shown). 
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Fig. 6.— Mass-age diagrams (MADs) for the components of V40 and V41. Results from the 
analysis of {M — R) and [M — L) diagrams are shown with solid and dotted lines, respectively 
(see Section [5] for an explanation). 














